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1 Introduction 

Mathematical models and numerical simulations for thin-film flows of a fluid 



have important applications in industrial and natural processes ]TT)| 



|13| , IfUfl , P, |J. Herein, we consider the slow motion of a thin liquid 
layer of an incompressible, Newtonian fluid over an arbitrary solid, stationary 
curved substrate. In a three dimensional and very slow lubrication flow, a 
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model for the evolution of a film on a curved substrate is shown II to be 



r] 2 (Vk - ]^i 4 (kI - K) • Vk 



where ( = r/ — ^nr/ 2 + |fciA>2?7 3 is proportional to the amount of fluid locally 
"above" the substrate; k is the mean curvature of the free surface of the film; 
K is the curvature tensor of the substrate; hi, k 2 , and k = k\ + k 2 are the 
principal curvatures and the mean curvature of the substrate respectively; 
and the operator V is defined in a coordinate system of the curved substrate 
as introduced in Section 2. This model accounts for the curvature of the 
substrate and that of the surface of the film. However, in many applications 
this model of slow flow of a thin fluid film has limited usefulness; instead 
a model expressed in terms of both the fluid layer thickness and the lateral 
velocity is needed to resolve faster wave- like dynamics 0, pi 10]. Roberts 
derives such a model for two dimensional flow. Here, based upon the Navier- 
Stokes equations for a viscous fluid, Section 3, we derive the following model 
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for three dimensional flow: 

^ = -V • fou) , (2) 

_(9u 7T 2 U 7T 2 .„ . TN U 

ot 4 7y 2 12 r/ 

where 7£ is a Reynolds number of the flow; g s is the component of gravity 
tangential to the substrate; and u = (mi,u 2 ) is the depth-averaged average 
lateral velocity. The first equation is a direct consequence of the conservation 
of fluid. The model's second equation incorporates inertia, viscous drag on 
the substrate, surface tension forcing caused by gradients of curvature, grav- 
itational forcing, and the geometric complexity of the substrate, respectively. 

The asymptotic accuracy of this model is assured by a systematic deriva- 
tion based upon centre manifold theory || as applied in Section 4. The 
physical fields associated with the above model are deduced as part of the 
procedure; approximations of the velocity and pressure fields are recorded 
in Section 5. The model @ is very general: it encompasses all substrate 
shapes and contains the lubrication model simply by setting TZ to zero and 
substituting for u in the first equation. This model describes the dynamics 
of a very wide range of thin fluid flows. 



2 The Orthogonal Curvilinear Coordinate Sys- 
tem 

Let S denote the solid substrate. If S has no umbilical point, i.e., there 
is no point on S at which two principal curvatures coincide, then there are 
exactly two mutually orthogonal principal directions in the tangent plane at 
every point in S Theorem 10-3]. Let ei and e 2 be the unit vectors in 
these principal directions, and e3 the unit normal to the substrate in the side 
of fluid flow. These basis unit vectors determine a curvilinear orthonormal 
coordinate system (xi, x 2 , y). Such a coordinate system is called a Darboux 
frame |J. The corresponding metric coefficients of the coordinate system are 

hi = mi(l - hy) , h 3 = l. 

Let y = r)(t,xx,X2) describe the free surface of the fluid. As derived by 
Roy et al [|11|, eqn(37)], the mean curvature of the free-surface is 
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where hi = mj(l— kiT)) are the metric coefficients evaluated on the free surface 
and 



A = sjh\hl + hl v l + hW X2 

is proportional to the free-surface area above a patch dxi x dxi of the sub- 
strate. 

We consider that the spatial derivatives of the fluid flow and the curva- 
tures of the substrate are much smaller than the metric coefficients of the 
coordinate system, so that an approximation to k, as needed for (|Ij), is 
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1 — kiT] 1 — k 2 T] 
where in the substrate coordinates 
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+ O (k 3 + V 3 r] 
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d ( m 2 drj \ d ( mj drj 



dx\ \midx\j dx 2 \m2dx2 



For later use, also observe that on the free-surface, unit tangent vectors t 1; 
t 2 and the unit normal vector n are 
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(^ei+^e 3 )/V^+< > 
-h 2 rj Xl ei - hiri X2 e 2 + hih 2 e 3 

(h 2 Vx 1 ) 2 + (hiVx 2 ) 2 + {hih 2 y 



3 Equations of Motion and Boundary Condi- 
tions 

Consider the Navier-Stokes equations for an incompressible fluid flow moving 
with velocity field u = (m, u 2 , v) and pressure field p. The flow dynamics are 
driven by pressure gradients along the substrate and caused by both surface 
tension forces, coefficient a, varying due to variations of the curvature of the 
free surface of the fluid, and a gravitational body force, g, of magnitude g in 
the direction of the unit vector g. Let the reference length be a characteris- 
tic thickness of the film H, the reference time fiH/a, the reference velocity 
U = a/fi, and the reference pressure a/H. Then the non-dimensional incom- 
pressible and Navier-Stokes equations are: 



V • u = 



K 



du 

~dt 



u • Vu 



-Vp + V 2 u + bg . 
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where 1Z = apH/p 2 is a Reynolds number characterising the importance of 
the inertial terms — it may be written as UH/ v for above reference velocity — 
and b = pgH 2 /a is a Bond number characterising the importance of the 
gravitational body force. 

We asymptotically solve these non-dimensional continuity and Navier- 
Stokes equations in the curvilinear coordinate system described in Section 2 
with the following boundary conditions 

1. The fluid does not slip along the stationary substrate, that is 

u = on y = . 

2. The fluid satisfies the kinematic free surface boundary condition 

dr] ui drj u 2 dr] 

~aT = v ~T7> T7T~ on y = r] . 

at hi ox 1 h 2 ox 2 

3. The stress across the free surface is caused by surface tension, in non- 
dimensional form 



— pn + r • n = nn on y — r) , 

where p is the fluid pressure relative to the assumed zero pressure of the 
medium above, and f is the fluid's deviatoric stress tensor evaluated 
at the free surface. 



4 A Centre Manifold Basis for the Model 

Using centre manifold techniques |2|], we derive a low-dimensional model 
of the dynamics described in detail by the non-dimensional continuity and 
Navier-Stokes equations in the general curvilinear coordinate system. We as- 
sume that the spatial extent of the fluid flow and the curving of the substrate 
occur on a much larger scale than the thickness of the fluid. Thus, e is intro- 
duced to parameterise the relatively small effect of these spatial variations 
and curvature, i.e., we introduce re-scaled * variables 

d d ^ 

To justify treating the gravitational forcing as small, we may introduce a 
parameter (3 such that the Bond number b = f3 2 ; then after adjoining % = 
to the dynamical equations centre manifold theory justifies an asymptotic 
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expansion in f3. As introduced by Roberts |7|, |9|, we also modify the tan- 
gential stress on free surface in order to establish three critical modes in the 
centre manifold model rather than the one natural critical mode of lubrica- 
tion theory. The modification is parameterised by 7; we construct a centre 
manifold for asymptotically small 7; then evaluating the results for 7 = 1 
recovers a model for the physical dynamics. Previous work [[7) has shown 
that such evaluation of these low-order expansions at 7 = 1 is accurate. 

For convenience we drop hereafter the "*" superscript on all re-scaled 
variables. The non-dimensional continuity and Navier-Stokes equations in 
the curvilinear coordinate system become the following governing equations 
(where i' — 3 — i): 
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dv ui dv u 2 dv dv , w? 

m + %^ + %dx- 2 +v dy- + emih -hi 
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+ ^ 



d ( hi dv 



dx? \ h? dx-, 



d (. , dv 
+ tt hih 2 — 
dy I oy 



dx\ \ hi dxi 



where the re-scaled scale factors are hi = mj(l — ekiy), and the boundary 
conditions become 

u = on y — , 

di] 

dt = 



Ui drj 
e- 



u 2 dr\ 
e— 



t, ■ r ■ n 



hi dxi h 2 dx 2 
mimim 2 Ui 
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ony = 7], 
on y = rj , 
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n ■ r ■ n = p + ft on y — r) . 



where 



I = V( e Mi*x) 2 + (eM* a ) 2 + (M2) 2 , 
and unit tangent vectors t 1; t 2 and the unit normal vector n are 

U = (hiei + er] Xi e 3 )/li, 

n = (-eft 2 j/ ai ei - eh^e^ + hji 2 e^)/l. 

Then by adjoining the trivial dynamical equations 

we get a new dynamical system in the variables u, 77, p, e, 7 and /3. The 
original system will be recovered by setting e = 1, 7 = 1, /3 = How- 
ever, the two systems are quite different from the view of centre manifold 
theory. Theory justifies treating all terms that are multiplied by the three 
introduced parameters as nonlinear perturbating effects in the new system. 

Then the linear part (in this new sense) of the governing equations and 
boundary conditions have solutions: Ui = v = p = 0,r] = constant; and 
independently v = p = 0, Ui oc sm^uy/rf) exp(At) where 

. u 2 

A = —-—7 , such that u = tanu; , 
Kr] 2 

Thus, the critical modes (with zero eigenvalue) are associated with varying 
thickness 77, and independently with the two shear modes U\ oc y. The sys- 
tem we consider has three critical modes and three trivial parameter modes, 
all other modes decay exponentially quickly. Also, the nonlinear terms are 
continuous at least. Therefore a low-dimensional model of the system is 
justifiably constructed by centre manifold theory. 

Denote the variables in the original system by v(i) = (77, ui, u 2 , v, p) . Cen- 
tre manifold theory guarantees that there exist functions V and G respec- 
tively describing the shape of the centre manifold and the evolution thereon, 
namely 

v(t) = V(r7,ui,w 2 ) , 



such that » 
dt 



V 

Ui 

u 2 



G(r), U!,u 2 ) 
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where dependence upon the constant parameters (e, 7,/?) is implicit in the 
above, and «j are depth-averaged velocities measuring the amplitude of the 
shear modes u,i oc y. The aim now is to find functions V and G such that 
v(£) are actual solutions of governing equations. We calculate V and G by an 
iteration using computer algebra ||. Suppose that an approximation V and 
G has been calculated and that we seek corrections V and G'. Substituting 

= G + G' 

into the governing equations then rearranging, dropping products of correc- 
tions, and using the linear approximation wherever terms multiply correc- 
tions, we obtain a system of equations for the corrections which is in the 
homological form 

CV' + AG' = R, 

where £ is the linear part of the governing equations and the boundary 
conditions, A is a matrix, and R is the residual of the governing equations 
using the reigning approximations, V and G. The procedure for solving 
the equations is as follows: first, choose G' such that R — AG' is in the 
domain of C; second, solve CV' = rhs with the given boundary conditions; 
then regard V + V and G + G' as the new approximations of V and G 
respectively. Repeat the procedure until the residual R becomes zero to the 
required order of error. Then the low- dimensional model has the same order 
of error by the Approximation Theorem in centre manifold theory. A 
computer algebra program (obtainable from the authors) is run to perform 
the computations. The key to the correctness of the results we report is the 
correct coding of the residuals within the body of the iteration. 

5 The Low-dimensional Model 

First, computing to low-order in the small parameters gives the following 
fields in terms of the parameters 

O / \ o 

p = —Ke — Vr/e — ijK^e +r](Y — l)g 3 e 
-tT'u • Vr? (§7 - I 7 y + 2 - 2Y) e 2 
-V • u(2 + 2Y - | 7 + | 7 F)e 2 
+0(e 3 + M 3 + /5 3 , 7 2 ) , 
u s = u (2F + | 7 r - 7 F 3 ) e 



v = V + V 
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such that — 

dt 
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U-2 



7 



+rj 2 (V^ + g s ) iy2 + i F 3 

480 ' ~ 240 ' 80 ' J 

-tjku (iy - r 2 + ±r 3 + 

+ V K ■ u (|y - Y 3 - - ^j7^ 3 

+^ 5 ) e 2 + 0(6 3 + « 3 + /3 3 )7 2 ), 

v = u • (F 2 + i 7 ^ 2 - |7^ 4 ) e 2 

-77V • u (F 2 + \lY 2 - \lY A ) e 2 

+0(e 3 + u 3 + (3 3 , 1 2 ) , 

where Y = y/rj, u = |u|, u s = (■u 1 ,-u 2 ) ) and O (e p + vfl + /3 m , 7 n ) is used to 
denote terms s which satisfy that s/ (e p + u q + (3 m ) is bounded as (e, u, (3) — >(), 
or s/ 7 n is bounded as 7 — >0. In the expression for p, the first and second 
terms are effects of the substrate curvature and free-surface, the third term is 
the correction for the first, the fourth term is hydrostatic, and the others are 
effects due to the motion of fluid. In the expression for u, the first line is the 
basic shear profile modified by boundary conditions, the second and third 
lines are modification due to forcing, and the others are effects of curvature. 
The expression for v expresses the vertical component of velocity is only 
dependent of the variations of free-surface and the other components. 
The corresponding evolution on this centre manifold is then 

^ = -e 2 V.(7 7 u)+0(6 3 + ^ + /3 3 , 7 2 ), 

^ = (V« + gs) (! + e 

+*T 1 «u (|7 " |) e 
+ V - l K .u(f 7 -3)e 
— 3r?~ 2 u7 

+0(e 3 + u 3 + (3 3 , 1 2 ) . 

The first equation gives the general expression representing mass conserva- 
tion for fluid considered. The first line in the right-hand side of the second 
equation is forcing, the second and third lines are curvature effects and the 
fourth line is drag. 

To recover a model of the original dynamics, we need to set 7 = 1. But 
as is apparent from the above, every coefficient in the centre manifold model 
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is a power series in 7. From earlier computation in [[7], we estimate that the 
radii of convergence of such 7 series are much greater than 1. Thus after 
computing to higher order in 7, we calculate every coefficient in the model 
from the first five terms in 7 series by setting 7 = 1. We also need to set e = 1 
which is valid provided the length-scales in the resolved dynamics are indeed 
much larger than the fluid thickness, that is, if the gradients are small enough. 
The low- dimensional model then becomes as given in the Introduction by (||) 
with errors O (V 3 + u 3 + (3 3 ). 

6 Conclusion 

This model conserves fluid and accurately accounts for the effects of the 
curvature of substrate, surface tension, gravitational forcing and fluid inertia. 

The low-dimensional model given in this paper is a simpler version. One 
may adjust a dynamical model to suit a particular application. 
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